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Abstract 

There exist lattice actions which give cut-off independent physical 
predictions even on coarse grained lattices. Rotation symmetry is re- 
stored, the spectrum becomes exact and, in addition, the classical equa- 
tions have scale invariant instanton solutions. This perfect action can 
be made short ranged. It can be determined by combining analytical 
calculations with numerical simulations on small lattices. We illustrate 
the method and the benefits on the d = 2 non-linear <r-model. 
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1 Introduction and Summary 



In many theoretical problems in physics, at a given stage, numerical analysis is 
necessary. Numerical problems require some regularization, most often different 
meshes are used. If the lattice is coarse, the numerical procedure is relatively easy, 
but the results are contaminated by the artefacts of the regularization. It is a difficult 
problem to remove these artefacts by making the meshes systematically finer and 
then extrapolating. 

We shall discuss the problem of artefacts in the context of lattice regularized 
asymptotically free quantum field theories. Removing the cut-off effects to reach 
the proper continuum limit is a central problem in this field [p]]. The method we 
describe, however, might find applications elsewhere, for example in solving partial 
differential equations. 

In the study of quantum field theories cut-off effects show up in different ways. If 
the lattice is coarse grained then all the results are distorted, independently whether 
short or long distances are involved. Similarly, even on a fine lattice, at distances of 
a few lattice units the lattice structure, rather than physics dominates the results. 

Removing the cut-off effects requires a careful analysis. Fig. 1 illustrates this 
process on the example of the running coupling in the d = 2 non-linear a-model @ . 
The renormalized coupling at the distance scale L is identified here with m(L)L, 
where m(L) is the mass gap in a periodic box of size L. For small L (L <C l/m^, 
where is the mass in the infinite volume) this quantity can be studied in per- 
turbation theory and shows the properties of an asymptotically free coupling: as L 
is changed, the coupling runs according to the perturbatively known beta-function. 
For larger values of L the scale dependence of the coupling, for example the relation 
between g(2L) = m(2L) ■ 2L and g(L) = m(L) ■ L, requires a numerical calcula- 
tion. A lattice is introduced with lattice unit a, and for a given g(L) (in Fig. 1 
g(L) = 1.0595) g(2L) is calculated on lattices with L/a = 5,6,7, ...,16 and the 
results are extrapolated to L/a = oo. This is a non-trivial numerical problem even 
in d = 2. 

A possibility to reduce lattice artefacts is to use improved actions. The idea 
behind this procedure is that for a large class of actions the physical predictions 
are universal while the cut-off effects are not. Although the history of improved 
actions is long ||, |J they had a limited effect on actual calculations until now. 
The reason is that in some cases the proposed actions were derived using ad hoc 
assumptions, or uncontrolled approximations, in other cases they were fixed by 
perturbation theory and it is not clear whether in typical simulations they improve 
and to what extent. Deep in the continuum limit a perturbatively fixed action 
certainly improves. This is, however, not necessarily the case at moderate correlation 
lengths. Fig. 1 illustrates this problem also. Perturbation theory fails to reproduce 
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even the sign of the observed cut-off effects |2[. 

A radical solution would be to use a perfect lattice action which is completely 
free of lattice artefacts. That such perfect actions exist follows from Wilson's renor- 
malization group (RG) theory [||, ||. The aim of this paper is to show that in 
asymptotically free theories a combination of analytical and numerical techniques 
allows finding the perfect action to a sufficient precision. We consider an action per- 
fect in practical calculations if even on a coarse grained lattice (correlation length 
—0(a)), no cut-off effects can be seen in numerical simulations. The perfect action 
is not unique. This fact can be used to find a perfect action whose range of interac- 
tion is short and whose structure is relatively simple. Simulating such actions the 
gain/cost ratio can be very large. 

In this work we study this problem in the d = 2, 0(3) non-linear cx-model, which 
we consider a pilot project, a preparation for d — 4 Yang-Mills theories and QCD. 
A quadratic lattice is considered and the partition function is written as 

Z = J] / dSJ(Sl - l)e^ (S) , (1) 

n J 

where (3A(S) is some lattice representation of the continuum action 

(3A cont (S) = ^J d 2 xd^S{x)d^S{x) . (2) 

Beyond the basic requirements of 0(3) symmetry, locality, correct classical limit, 
translation and 90°-rotation symmetry, the form of the lattice action is largely 
arbitrary. It might contain nearest neighbour, next-to-nearest neighbour, etc., 
even different multispin interactions. Let us denote the corresponding couplings 
by Ci, c 2 , . . .. The action /L4(S) is represented by a point in the infinite dimensional 
space of couplings (/3, ci, C2, . . .). We shall consider RG transformations in configu- 
ration space, namely block transformations with a scale factor of 2. Under repeated 
block transformations the action moves in this coupling constant space. The ex- 
pected flow diagram is sketched in Fig. 2 || |7] . In the (3 = oo hyperplane there is a 
fixed point (FP) c*, dj, . . ., whose exact position depends on the details of the block 
transformation. We shall use the notation *4(S; c\, c\, . . .) — A*(S) and call /L4*(S) 
the FP-action. The FP has one marginal and infinitely many irrelevant directions. 
The marginal operator is A* itself [[5J. Actually, A* is not exactly marginal, it is 
weakly relevant. If (3 is very large, under a RG transformation with a scale factor 2, 
(3A*{S) — > /3'A*(S), where (3' = (3 — In 2. The coupling g = 1/(3 grows according 
to the beta-function of this asymptotically free theory. The trajectory which leaves 
the FP along the weakly relevant direction is called the renormalized trajectory 
(RT). For large (3 the RT runs along the FP-action, but for smaller (3 they do not 
coincide anymore. 

It is easy to see that the points of the RT define perfect actions. The argument 
goes as follows. At any given f3, the point of the RT is connected to the infinitesimal 
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Figure 2: Flow of the couplings under RG transformation in the O(N) non-linear 
a-model 

neighbourhood of the FP by (infinitely many steps of) exact RG transformations. 
Since each step increases the lattice unit by a factor of 2, any distance at the given f3 
(even 1 lattice unit) corresponds to a long distance close to the FP. The infinitesimal 
neighbourhood of the FP is in the continuum limit, there are no cut-off effects there 
at long distances. On the other hand, for all the questions which can be formulated 
in terms of the degrees of freedom after the transformation we get the same answer 
as before it. So, there are no lattice artefacts at the given (3 on the RT at any 
distances. 

Since a RG step is a non-critical problem, the FP-action and the actions on 
the RT in general are expected to be local Locality, however, allows non- 

negligible interactions over several lattice units in the action. For practical reasons 
we neeed more than locality. For practical applicability we have to answer positively 
the following questions: i. Is it possible to determine *4*(S) to a good precision? ii. 
Is «A*(S) of sufficiently short ranged? Is the structure of A* simple enough allowing 
a parametrization where the number of couplings remains relatively low, 0(10-100)? 
iii. Questions i. and ii. for the points of the RT. 

Our pilot study shows that for the non-linear a-model, using a properly chosen 
RG transformation, the answer to these questions is 'y es \ even if we go down to 
small correlation lengths . The determination of ,A*(S) is a saddle point problem 
which requires minimalization over classical fields. The FP, ^4*(S) is the perfect 
classical theory on the lattice: it has stable instanton solutions whose action is 
independent of the size p of the instanton (scale invariance). This should be the 
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same in d — 4 gauge theories also. This fact might open new ways in the study 
of topological effects on the lattice. The range of interaction in *4*(S) depends 
on the RG transformation. By a proper choice of the block-transformation *4*(S) 
becomes surprisingly short ranged. ^4*(S) contains multispin couplings also but its 
structure is relatively simple. With O(20) couplings an excellent parametrization 
can be obtained which works even on coarse configurations, i.e. at small correlation 
lengths. The problems related to *4.*(S) can be solved partly analytically, which is 
a special bonus in asymptotically free theories. 

If j3 is not very large, the FP-action Afp = /3fpA*(S) is not a perfect action. 
It is, however, the perfect classical action and we expect that the lattice artefacts 
are significantly reduced. Actually, we were not able to identify any cut-off effects 
when simulating Afp in our pilot study. The lattice artefacts in the problem of 
the running coupling in Fig. 1 can be used, for example, to illustrate how well the 
FP-action performs. As Fig. 3 shows, even the coarsest discretization L/a = 5 gives 
a result which is consistent with the final extrapolated value. Let us emphasize at 
this point that Afp has nothing to do with perturbation theory. 

The FP-action is 'almost perfect', which suggests that the RT runs close to 
the FP-action even at moderate /3-values. In this case the blocked FP-action will 
converge rapidly towards the RT. We illustrate on an explicit numerical experiment 
that the problem to find and parametrize the action after a block-transformation 
at finite (3 is feasible. In traditional Monte Carlo simulations one goes deep in 
the continuum limit, simulates large lattices and fights with memory problems and 
critical slowing down. Our message is: do careful calculations on small lattices 
instead. 

There remained many interesting problems even in this d = 2 model which we did 
not touch: an explicit study of the classical solutions, the possibility of using strong 
coupling expansion rather than Monte Carlo simulations for the perfect action at 
small correlation lengths and, most notably, the problem of constructing currents 
corresponding to the perfect action. 

The idea of constructing improved actions using RG ideas is old ||, [J. But, 
perhaps because the problem was considered as not really feasible, already at the 
start crude approximations and assumptions were introduced. We would like to 
show here that — due to the specific properties of asymptotically free theories - 
one can get close to the dream of a perfect regularization. First, finding the FP is a 
relatively easy classical problem. This is a significant advantage over other statistical 
models, e.g. the d = 2 Ising model |24|] . Second, even a simple optimalization makes 



the FP-action short ranged. Third, the FP-action, which is the perfect classical 
action, performs very well even at small correlation lengths. This feature makes it 
a good first approximation in finding the RT. 

Earlier RG and Monte Carlo RG studies M, M revealed many analogies between 
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d — 4 Yang-Mills theories and d = 2 asymptotically free spin models. This gives a 
hope that the results of this paper can be generalized. Of course, the final goal is to 
treat full QCD where the cut-off effects are even more troublesome. The formalism 



can be extended to fermions. The free fermion problem has been discussed in fLO 
but the real hard work is to be done. 



We should mention that the physical picture in the multigrid-RG program as 
formulated by Mack [23] has close analogies to that discussed here. 



2 The fixed point and the fixed point action 



2.1 The equation for A* 



We shall define the RG transformation in the model in eq. ([!]) as follows. The square 
lattice is divided into 2x2 blocks and to every block we associate a block spin R nB , 
which is a certain average of the four original spins in the block. We shall denote 
the points of the original and the block lattice by n and ris, respectively. The spins 
S n and R n£j are normalized 0(N) vectors. The RG transformation is defined as 



-I3'A'(R) 



JS { n B 



PR nB ]T S n -lnY N (P\ Sn|) 

nans ridnB 



(3) 



/ = n / dS»* (S* - 1 



where P is a parameter of the transformation, while the normalization factor Yjv(z) 
assures that the partition function remains unchanged up to an irrelevant constant: 



dR5 (R 2 



1 e 



Rb 



const • Yi 



N 



(4) 



The function Y N (z) is related to the modified Bessel function (some of its prop- 
erties are summarized in Appendix F in [11|, for example), specifically Y%(z) ~ 
sinh(z)/z. For P — ► oo, the block-transformation goes over to a 5-function giving 



R 



J2nen B Sn/| J2n€n B S n |- We keep, however, P finite and write 

P = 0.[« + 0(1/0)] , 



(5) 



where k is a free parameter. Similar RG transformations were considered before in 
Monte Carlo RG studies [0. The coupling (3 is defined by the requirement that 
on very smooth configurations the action goes over to the form in eq. (Q) (see also 
section 2.4). 
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In the limit (3 — > oo, eq. 

e -^'(R) = exp L (s) _ R £ 



can be written as 
R 



»B 



■"3 E ^ n I E ^' 



(6) 



where /3' = /? — 0(1), due to asymptotic freedom. Eq. (|) is a saddle point problem 
in this limit, giving 



^(R)=nunU(S)-«X; 



»B 



n£n B n£n B 



(7) 



The FP of the transformation is determined by the equation 
A*(R) = mm <A*{S) 



«B 



^™B ^ n I ^ r 



As we discussed at length in the Introduction, the FP-action /L4*(R) - although it 
is not 'perfect' at finite (3 - is an excellent starting action for further studies even 
at moderate /^-values (at moderate correlation lengths). So, we shall need *4*(R) 
for strongly fluctuating configurations also. Observe that eq. @ reduces to the 
saddle point equation eq. (0) for any configuration R. If the configuration R is 
strongly fluctuating, then the minimizing configuration S will not be smooth either. 
In general, eqs. (0),@ and their solutions have nothing to do with perturbation 
theory. 

Eq. (|8D for A* is a highly non-trivial equation. Some of the important properties 
of the solution can be obtained, however, without solving the equation explicitly. 



2.2 A* as a perfect classical theory on the lattice 

Using eq. (§) it is easy to show the following^ 
Statement 

If the configuration {R} satisfies the FP classical equations (i.e. the classical equa- 
tions corresponding to A*) and it is a local minimum of A*(R) (allowing also for 
zero modes) then the configuration {S(R)} on the finer lattice which minimizes the 
right hand side of eq. (|8|) satisfies the FP equations as well. In addition, the value 
of the action remains unchanged, ^4*(S(R)) = *4*(R). 

Proof 

Since {R} is a solution of the classical equations of motion, 5A*/5H = 0, the 
configuration {S} = {S(R)} should satisfy the equation 



E Sn 

n£n B 



(9) 



3 We are indebted to Uwe Wiese who raised our attention to the classical solutions in this 
context. 
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for any Here A„ fl > since we excluded the case when *4*(R) is a local maximum 
Since the expression 



in R nfl 



II B 



nGris nGn B 



(10) 



in eq. (§) takes its absolute minimum (zero) on the configuration {S(R)} satisfying 
eq. (|D|), {S(R)} is also a stationary point of A*, 

5A*(S) 



= o, (n; 

S=S(R) 



5S 

with the same value of the action, ,A*(S(R)) = A*(R), what we wanted to show. 



According to this Statement, to any solution {R} (which is a local minimum) 
with a characteristic size p, there exists another solution {S(R)} of size 2p with 
exactly the same value of the action A*. It is natural to assume that this solution 
on the finer lattice is also a local minimum. Specifically, for N = 3 this observation 
implies that if A* has an instanton solution of size p then there exist instanton 
solutions of size 2p, 2 2 p, . . . , 2 k p, . . . with the value of the action being exactly Air for 
all these instantons. The value 4n follows from the fact that very large instantons 
are smooth on the lattice and then any valid lattice representation of eq. (0) gives 
the continuum value. 

The standard action has classical solutions of antiferromagnetic character — e.g. 
when one of the spins points downwards in a background of upward pointing spins 
- which are smooth local maxima in some of the spins. We expect, although we 
are not able to prove, that A* does not have such solutions. For example, the 
configuration mentioned above is not a solution since 5A*(R)/5H does not vanish 
at •& = 7T, as shown in Fig. 6. 

It was observed earlier in [jnj that using an approximate RG improved action 
the stability of instantons can be increased significantly. This result is in accordance 
with the Statement above. 

It is important to observe that the reverse of the Statement is not true: if 
the configuration {S n } is a solution, then the configuration {R ns }, where R rifl = 
J2nen B S n /| J2nen B S n |j is n °t necessarily a solution. The proof fails because for 
this configuration {R„ B } the minimizing configuration {S(R)} in eq. (0) is not 
necessarily equal to {S n } itself. One can show that {S n } is a minimum but not 
always the absolute minimum. Actually, this is the mechanism which prevents the 



existence of arbitrarily small instanton solutions on the lattice finfl 



With respect to the classical solutions, A* has the same scale symmetry as the 
continuum action. We shall give later further arguments which show that it is a 
perfect classical lattice representation of the continuum action. 
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2.3 Parametrization 



Eq. (|D defines the value of the FP, A* for any configuration R and, as we shall 
discuss later, one can advise numerical procedures which give this value to high 
precision. However, if we want to use the FP action (3ppA* in numerical simulations 
we can not avoid to parametrize A* by a (limited) number of coupling constants. 
We can write in general 

A*(S) = -^5>(r)(l-S n S n+r ) (12) 

* n,r 

+ J2 c(7n,n2,n 3 ,n 4 )(l-S niS n2 ) (1 S raa S n4 ) -\~ . . . 

711,712,713,714 

where the summations go over all the lattice points. It is a significant help in the 
parametrization and optimalization problem that the first two functions p and c in 
eq. (|T2| ) can be calculated directly. We shall discuss this problem first. 



2.4 The determination of p(r) 

Consider a configuration R where the spins fluctuate around the first axes: 

Rn B = N 1 '^* | (13) 
V Xn B 

where Xn B has (N — 1) components, and \Xn B \ <C 1- In this case, the saddle-point 
configuration in eq. (ffi) will also fluctuate around the first axes, so we can write 



1 - 7T 2 

n 



7TV 



(14) 

where |7? n | <C 1. The terms quadratic in tt and \ gi ye a closed equation for p when 



eqs. (PD and (O) are used: 



9 Z! P( r B)Xn B Xn B +r B = mjn J ^ p^T^TT^ + 2k £ { Xn s ~ T 2 J [ ■ 
Z n B ,r B W [ Z n,r n s \ 4 n£n B / J 

(15) 

The equation for p is independent of N as far as N > 1. Eq. (|l|) is just the FP 
equation of a free scalar theory with a Gaussian block-transformation: 



C " eX p|-^ Z P( r B)Xn B Xn B +r B \ = 
{ ^ n B ,r B ) 

II / d ^ exp < -2K Z ( Xn B ~ J ^» ) ~ ^ Z PM^n+r | 

n I n B \ 4 n£n B / z n,r J 



(16) 



Eqs. ( p~5| ) and ([X6|) are equivalent: Gaussian integrals can be replaced by minimal- 
izations. The free field problem has been studied and solved in |Tj| a long time ago. 
In Fourier space the solution has the form 

_L = y 1 tt si " 2(9 " /2) + J- (it) 

where the summation is over the integer vector I = (/ , h), (q + 2nl) 2 = (g + 27r/ ) 2 + 
(qi + 27r/i) 2 and 

(18) 

The normalization of p is fixed by demanding p(q) — > q 2 for small q. In config- 
uration space this corresponds to J2 r p(r)r 2 = —4. We fix the value of p(0) by the 
convention J2 r p( r ) = 0. 



There is an elegant, simple way to derive eq. (|TT| ) which we discuss below, since 
it might simplify the analogous problem for gauge and fermion fields also. A similar 
trick was used in |TJ| in an approximate RG calculation in the large-iV limit with 



^-function blocking. 

The observation is that in a free field theory the two-point function is directly 
related to the action. On the other hand, the two-point function at the FP can be 



calculated easily. From eq. (16) it follows 



/1\ 2 1 

(Xn B Xn' B ) = ( t) E + J^n B n' B ■ (19) 



n'Gn' 



Using this relation recursively, after j RG steps we get 

(Xn B Xn> B ) = ( j) ' E i^n') + ^( 1 + i + ^ + -- - + i7) 6 nsn' B ■ (20) 



n£n B 
n'en' 



where x is the block-field after j steps, and the block hb contains 4 J original lattice 
points. Taking the j — > oo limit, and considering the original lattice infinitely 
fine, the summations in the first term on the right-hand side of eq. (|20D go over to 
integrals. The integrals are over a block (a x a) centered around the points hb and 
n' B , respectively. For the two-point function on the infinitely fine lattice we can use 
the standard 1/k 2 propagator. Taking a = 1, and simplifying the notations we get 

rl/2 rl/2 roo (pk Q ik ( n + x - n '- x ') I 
= Ll/2 ^ 1-1/2 ^ loo J2^ 2 k 2 + ■ (21) 

After many RG steps the system runs to the FP (taking the mass to be zero at the 
beginning, only irrelevant directions remain), therefore 

r d 2 q e 9(n - n,) . , 

b^ )= Liw~wr- (22) 
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Dividing the /c-integration into an integration and summation by introducing 
the integer vector I = (/ , li) 



A.; 



2irl h 



q.i e (— 7T,7T) 



0,1 



(23) 



and performing the integrals over x and x', the result in eq. (|17]) follows. Using 
similar steps, the two-point function in configuration space has been obtained for 
^-function blocking by Iwasaki ][| and by Gawedzki and Kupiainen [17] also. 



Let us remark that the 1/4 factor in the block transformation of a free theory 
in eq. ([16]) is replaced by b = 2~ (c(+2 )/ 2 in d dimensions. Only this choice leads to a 
FP || 15]. In the non-linear cx-model, the factor 1/4 in the last term in eq. ([T5|) is 
fixed. That these two factors coincide for d = 2 is related to the fact that in two 
dimensions (and only in this case) the non-linear cx-model is asymptotically free. 



2.5 The properties of p(r); the perfect Laplace operator on 
the lattice 

The properties (most notably the range of interaction) of p(r) depend on the free 
parameter k. The choice k ph 2. is optimal. For k = 2, p(r) decays very rapidly with 
growing |r|, like ~ exp(— 3.44|r|), Fig. 4. For this choice p(r) is strongly dominated 
by the nearest neighbour and diagonal couplings while the couplings at distance > 1 
are already small. As Table 1 shows, p(3, 3) is, for example, 5 orders of magnitude 
smaller than the nearest neighbour coupling. The k = 2 value is distinguished also 
by the fact that in d — 1 (or, equivalently, in any d on configurations which depend 
only on one coordinate) this choice leads to a simple nearest neighbour action. Fig. 4 
also shows that the choice k = oo (corresponding to a block transformation with 
5-function) gives a considerably larger interaction range, p(r) ~ exp(— 1.45|r|). 

A free field theory defined in terms of p on the lattice gives a cut-off independent 
account of the physical properties of the theory. The spectrum is exact. In the 
continuum the energy of an excitation with momentum k\ is E{k\) = \k\\, where 
k\ e (— oo, oo) (massless relativistic free particle). Using eq. flT7| ) one can check easily 
that the lattice field theory with p has exactly this spectrum. (A possible method 
is to calculate the two-point function at large time-separation.) Although qi in 
eq. ( |17D lies in (— n, n) only, the full spectrum is reproduced because for every given 
qi there are infinitely many energy eigenstates. This is related to the summation 
over / in eq. (|17l). 

Fig. 5 compares the predictions for the energy spectrum of p (exact) to those 
of the standard lattice regularization, Symanzik tree level improved action and the 
approximation where the p(l,0) and p(l, 1) couplings are kept only.Q Already this 
crudest approximation performs very well. 

4 We are indebted to M. Blatter for his help on Fig. 5. 
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r 


p(r) 


r 


p(r) 


(1,0) 


-0.61802 


(4,0) 


-2.632 • 10" 6 


(1,1) 


-0.19033 


(4,1) 


7.064 • 10- 7 


(2,0) 


-1.998- 10" 3 


(4,2) 


1.327- 10^ 6 


(2,1) 


-6.793 • 10~ 4 


(4,3) 


-7.953- 10" 7 


(2,2) 


1.625 • 10~ 3 


(4,4) 


6.895- 10- 8 


(3,0) 


-1.173 ■ 10~ 4 


(5,0) 


-8.831 ■ 10~ 8 


(3,1) 


1.942 • 10~ 5 


(5,1) 


3.457- 10- 8 


(3,2) 


5.232 ■ 10~ 5 


(5,2) 


3.491 ■ lO" 8 


(3,3) 


-1.226 ■ 10~ 5 


(5,3) 


-3.349 ■ 10~ 8 






(5,4) 


8.408 • 10" 9 






(5,5) 


-1.657- lO" 10 



Table 1: The couplings of the quadratic two-spin interaction terms at a distance 
r = (r ,ri) for the optimal choice of the block transformation with k = 2. Note 
that in our convention, for the standard action the only non-vanishing entry in this 
list would be psrCLO) = —1. 



The finite difference operator Ai&) n = — J2r p( r )$n+r is a perfect discretization 
of the continuum Laplacian. Let <p(x) be a field in the continuum satisfying the 
Laplace equation 

A(j>(x) = . (24) 

Lay a lattice on the d = 2 continuum plane with lattice unit 1 and define on the 
lattice point n a lattice field 

$ n = d 2 x<f)(n + x) (25) 

"2 

The field $ n satisfies 

A L $ = , (26) 

and gives the same action as 4>(x): 

- J rf 2 x0A0 = -$>„ (A L $) n (27) 



2.6 The determination of (2(711,712,713,714) 

In order to obtain the coefficients c(. . .) we have to expand eqs. (|],[12]) up to quartic 
order in 7? and \. In quadratic order, discussed before, the minimizing it configura- 
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tion in eq. ([§]) is linear in \'- 



n n = Y<M(n,n' B )xn' B , (28) 
n' B 

where the matrix M is short ranged (if n is in the block n B then only small \n B — n' B \ 
gives significant contribution in the sum in eq. Q2"8"|)) and can be obtained explicitly 
when solving the free field problem ]T3]. In order to obtain the \ dependence of 



the right-hand side of eq. (§) up to quartic order it is sufficient to substitute the 
minimizing value of ir up to linear order, eq. (|28|) . The cubic corrections to eq. ([28" ) 
give 0(x 6 ) terms only in eq. (|8|). 

In quartic order eq. (|8]) can be written in the form 

H A ( n B, n' B ; n" B , n" B ) {xn B Xn' B ) (xn B Xn%) = , (29) 



where the matrix A depends linearly on c(. . .) and is symmetric under % «-> n' B , 
n" B n' B , and (riB,n B ) <-> (n B ,n B ). Its form is not very illuminating and can be 
easily derived, so we do not quote it here explicitly. Although the matrix A has no 
explicit iV-dependence, the solution c(. . .) distinguishes between N > 3 and N = 2. 
For iV > 3 eq. (E^) implies A = 0, while for N = 2 \ is a 1 -component vector and 



eq. (|29|) gives 

A(n B , n' B - n" B , n%) + A{n" B , n' B ; n Bl n%) = , (N = 2) . (30) 



These equations can be solved for c(. . 
dure. The short range nature of p(r) for 
Table 2 enumerates those couplings with 



.) by a rapidly converging iterative proce- 
k = 2 is inherited by c(. . .), as expected. 
\c{...)\ > 5-1Q- 5 . 



2.7 The numerical procedure to solve the FP equation 

For a given configuration {R}, eqs. (0,§) can be solved numerically for the FP value 
^4*(R). Calculating *4*(R) for a large number of different {R} configurations, the 
important couplings (beyond those in p(r) and c(. . .)) can be found by a fitting 
procedure. The finite size distortion on the couplings will be small if the lattice on 
which {R} is defined is larger than the range of interaction in A*. Since the later is 
small for k = 2 (see p(r) and c(. . .), and further discussion in 2.8) this requirement 
is easy to satisfy. In this pilot study we worked on 5 x 5 periodic lattices when 
looking for the couplings of the FP, A* numerically. 

We worked out two different numerical procedures to calculate ^4*(R) for a given 
configuration {R}. Both were based on iterating eq. (^|). For a process with k 
iterations we define the lattice sizes L (0) , L w = L^/2,..., = L { ^/2 k . The 
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# type 



coupling 



# type 



coupling 



11 



13 



J 

X 




JL 




0.05344 
0.02155 
0.01209 
0.00144 
0.00121 
0.00062 
-0.00057 



10 
12 
14 




0.00960 



0.01881 



6 * * -0.00258 




I / 



0.00123 



-0.00070 



0.00062 



-0.00042 



Table 2: The coefficients of the leading quartic terms (1 — S ni S n2 )(l — S„ 3 S n4 ) of A*. 
A term which is defined by the sites ni, ri2, 713, 714 and the topology of the four-spin 
coupling occurs only once with these coefficients in the action. The numbers in this 
Table are multiples of c(nj, . . . ,n 4 ) in eq. (|T2|). This symmetry factor is 4 for the 
couplings #1,2 and 8 for the others. 
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configuration {R} lives on the lattice of size L( fc ) (=5, in our case), while we denote 
the spins on the th , 1 st , . . . , (k - l) th level by S<g, S«,. . . , Sg^J. We get 



^( fc )(R)= min \A {0) (sW)-kY 
{s( ),s( 1 ),...,s(*-i)} \ y 1 

712 



c(i) v S(°) - I V S(°)| 



noGni 



no6m 



e(2) V - I V 
°n 2 Z^ °ni I Z^ °m 



nie?i2 



riiGn 2 



(31) 



E 



E s?-? -I E sit 



e-l) 

'rifc_i i "n k -i 
n k _ 1 £n k n k _ 1 £n k 



Here .4.( ) is the zeroth order approximation for A*, while A^ k '(R) is the k th approx- 
imation for the number *4*(R). For A^ ' one might take even the standard action, 
but the convergence in k will improve if a better approximation is used. (For exam- 
ple, one can use the results on p(r) and c(. . .), or those of earlier numerical runs.) 
In eq. ( [5ID Dn eni denotes a sum over the 4 spins in the block m, etc. 



Note that eq. (pl|) and the way we solved it are strongly reminiscent to a multigrid 



procedure [^3], |T 



In the first method the minimalization in S was performed locally and the lattices 
were swept through until A^ (R) became stable to a given precision. In the second 
method we used annealed cooling. The largest irrelevant eigenvalue of the FP is 1/4 
in agreement with the observed convergence rate in the number of blocking steps k. 
Even taking the crude approximation for A^ where only the nearest neighbour and 
diagonal couplings from p(r) are kept, |^4^(R) — ^4*(R)|/|^4*(R)| was estimated 
less than O(10 -4 ). The two different methods gave consistent results. 



2.8 The properties of A* on coarse configurations 

On smooth configurations the first terms in eq. QT2j) dominate. These terms are 
short ranged and explicitly known. It is important to investigate, however, whether 
A* remains short ranged even on coarse configurations. We need further hints on 
the general parametrization problem also. 

In order to see how *4*(S) behaves as the relative angle between the spins is 
increased, we rotated a single spin S away from the I s * axis by an angle while all 
the other spins pointed in the I s * direction. The corresponding value of A*, which 
we denote by Ai(&), is shown in Fig. 6. The curve is very simple and can be fitted 
easily if we parametrize it in powers of d 2 /2 rather than in powers of (1 — cost?). 
This empirical observation is also supported by considering a quasi one-dimensional 
solution to the equations of motion, with S(n) = (cos^rio, sin^no, 0), where the d 
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Ai(#) 




Ai{$) - \A 2 {-&) 


-p(2,0)tf 2 /2 


0.09967 


0.016094 


0.032166 


0.000011 


0.000010 


tt/4 


0.99558 


1.98986 


0.00065 


0.00062 


tt/2 


3.9339 


7.8627 


0.0025 


0.0025 


2.678 


10.9999 


21.9863 


0.0067 


0.0072 


2.85 


12.3352 


24.6553 


0.0075 


0.0081 



Table 3: One and two spins (at a separation (2,0)) are rotated in a trivial back- 
ground. The values of A* for the corresponding configurations are given together 
with Aifi) - \A 2 (■&) which measures the direct (2,0) interaction. The last column 
shows the contribution from the first term in eq. (12'). 

dependence of A* is exactly given by 2 /2. So, we rewrite eq. ( |T2"D as 
-4( S ) = -^Ep( r )^<„+r+ E c(n 1 ,n 2 ,n 3 ,n 4 )^< in2 -^<3 in4 + ... (12') 

n,r z ni,rt2,n,3,n4 ^ ^ 

where $ niin2 is the angle between the spins S ni and S„ 2 . The function c is trivially 
related to p and c of eq. ([12]) . Using this parametrization, the leading couplings of p 
and c will dominate A* even on coarse configurations. This is illustrated in Fig. 6, 
where a did 2 /2-\-d2{'d 2 /2) 2 form, with d\ and g?2 obtained from p and c in Tables 1,2, 
agree very well with the data. 

Next we wanted to check whether the distant couplings remain small even at 
large relative angles. In a trivial background (spins in the 1 st direction) we rotated 
2 spins, which were at a distance A from each other, by an angle $ relative to the 
I s * axes (keeping them parallel to each other). Denoting the corresponding value of 
A* by A 2 ($) it is easy to see that A 1 ($) — |A 2 ($) is a measure of the direct coupling 
between the two spins. Table 3 gives the measured values at different $ values for 
distance A = (2,0) compared with — p(2, 0) $ 2 /2, the first, exactly known term in 
eq. (12'). The direct coupling remains small and very well given by this single term 
even at large angles. 

2.9 Fixing the couplings in this pilot study in 0(3) 

In the following we consider the 0(3) model, N = 3. We generated ~ 500 different 
R configurations and calculated the corresponding values of the FP, ^4*(R) using 
the minimalization procedure discussed in 2.7. Among these configurations there 
were smooth, but also quite coarse ones, typical for a correlation length of a few 
lattice units only. We made fits with up to 69 different couplings, taking them all 
short ranged, with the largest distance being (2,2). 
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# type coupling 
1 • • 0.61884 

0.19058 

0.01881 



type 



10 
13 
16 
19 

22 



A 

X 



0.02155 
0.01078 
0.01209 
-0.00258 
-0.01817 





coupling 
-0.04957 

-0.02212 

-0.00139 
0.00717 
0.00765 

-0.00114 
0.00387 

-0.00772 



type 







coupling 
-0.01163 

-0.00463 
0.00497 

-0.00055 

-0.00557 
0.00548 

-0.00100 
0.04970 



Table 4: The couplings used for the FP-action in eq. (12'). The notation is explained 
in the text. The coefficients of the quadratic and quartic interactions are calculated 
analytically, the higher order interactions are fitted. 



Fig. 7 shows a fit with 24 couplings, 8 fixed analytically by p and c, 16 fitted. 
The couplings and their fitted values are given in Table 4. The couplings #1 and 4 
are slightly renormalized to satisfy the condition Y^, r p{r)r 2 = —4 exactly. In this fit 
all the couplings can be put on a 1 x 1 square. In Table 4 the couplings are given 
in a graphical notation. The coupling #15, for example, multiplies the expression 
~$ 2 A ■ ~$ 2 B ■ (\'&c) 2 1 where = acos(SiS2), $b = acos(SiS3) and $c = acos^Ss). 
Here Si, S2, S3 are the 3 spins forming the triangle, Si sitting at the 90° corner. 

As Fig. 7 shows the fit becomes poorer on coarse configurations (those have 
larger action values), in a few cases the error reaches ~ 1%, but for most of the 
configurations the error is much smaller. Including the (2,0), (2, 1) and (2,2) cou- 
plings the fit becomes somewhat better, but for this pilot study we decided to use 
the 24-parameter fit in Table 4. 

We want to emphasize two points here. Using 24, or even twice as many couplings 
in a Monte Carlo simulation is not really a problem - at least in the spin model 
considered here. Although this action includes 3- and 4-spin interactions also, one 
can generalize the cluster Monte Carlo technique |0| for this case. Our cluster 



16 



program runs ~ 3 times slower with 24 parameters than with the standard action. 
Relative to the gain we are looking for, this is irrelevant. The second remark concerns 
the numbers in Table 4. The fitted parameters are effective values, one must not 
associate significance to a single entry in this Table. Modifying somewhat the set of 
operators used in the fit, these entries would change, but the global quality of the 
fit would hardly be different. 

2.10 Results obtained by simulating the FP action 

The FP-action A* FP is defined as (3fpA*(S). The coupling j3 FP is fixed uniquely by 
the normalization condition on p (sect. 2.4). For large @fp, the RT runs together 
with the FP-action (Fig. 2), therefore A* F p is the perfect action. For intermediate 
(5fp values it is not perfect anymore, but - as we shall see - it performs amazingly 
well. As we discussed in the previous sections, the FP-action can be determined 
with the help of classical calculations. For n = 2 it is very short ranged and a 
relatively simple parametrization describes it well. This parametrized form (given in 
Table 4) can be easily simulated. Our cluster Monte Carlo program performed very 
effectively. We would like to urge workers in this field to replace the standard action 
by the FP-action as a first approximation towards 'perfectness'. These remarks 
apply to the CP n -model also, especially in studies related to topology (sect. 2.2). 

In the following, for simplicity, we shall call the 24-parameter action in Ta- 
ble 4 the FP-action keeping in mind possible effects from the errors of this simple 
parametrization. In the Appendix we calculated the relation of the A-parameter of 
the FP-action to that of the standard action. The result 

— — = 8.17, (N=3, couplings from Table 4) . (32) 

A S T 

can also be written as 

Pfp = Pst - 0.334 + 0(1//? 5 t) • (33) 

Eq. (|33D gives a feeling about the meaning of /3fp- At intermediate (3 values we 
give — for orientation — the infinite volume correlation length £,(/3fp) at a few (3pp 
values: f (1.214) « 36, f (1.08) « 18, f (0.85) « 6. 

2.10.1 Cut— off dependence of the running coupling 

We discussed this problem in the Introduction already. We have chosen the g(L) = 
m(L)L = 1.0595 point from ref. since here the observed cut-off effects in the 
relation m(2L)2L versus m(L)L were large and opposite in sign with respect to 
perturbation theory (Fig. 1). We simulated the FP-action on a lattice of spatial 
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size L/a = 5 and tuned (3fp until m(L)L became close to the prescribed value: 
at Pfp — 1-0821 we got g(L) = 1.0578(5). Then measuring on a lattice with 
L'/a = 2L/a = 10 at the same (3pp value we had m(2L)2L = 1.2611(9). In the 
time direction the lattice was chosen to be at least 6-times larger than the finite 
box correlation length = l/m(L) and distances larger than three times £ were 
used in the fitting procedure to obtain the mass gap. At the end we shifted g{2L) 
according to the slight difference between the actual (1.0578) and prescribed (1.0595) 
value for g(L) obtaining g(2L) = 1.2638(12). For this shift we used the universal 
curve obtained in ||. The error from this procedure is negligible. We repeated this 
calculation at L = 10 also, for (3pp = 1.214 with the results: m(L)L = 1.0613(8) 
and m(2L)2L = 1.2664(18). After shifting this gives g(2L) = 1.2635(22). In Fig. 3 
these two points are compared with the extrapolated prediction g(2L) = 1.2641(20). 
No cut-off effects can be seen.Q 

2.10.2 The two— point function 

We want to compare the spin-spin correlation function of the standard nearest 
neighbour model with that of the FP-action on coarse lattices and check the extent 
of violation of rotation symmetry. In this respect one has to understand the fol- 
lowing point. The correlation functions do not describe directly physics since they 
depend on the way the fields are defined. This is reflected, in general, in the scheme 
dependence of the wave-function renormalization. Here the relation between the 
fields of the FP-action and the perfect fields of the continuum is more involved. In 
the case of a free field, as the derivation in section 2.4 (especially eq. (pTD) clearly 
shows, the field of the FP-action is the integral of the perfect field over a square 
of length a. This definition brings an inherent rotation symmetry breaking in the 
two-point function. Similarly, in the full 0(3) model the shape of the block used 
in the RG transformation gives a non-rotation invariant definition for the fields in 
the FP-action. Of course, physical predictions, like the energy values related to the 
exponential decay of the two-point function, are rotation symmetric. 

The effect discussed above is, however, very small. This can be checked analyti- 
cally in the free field case and it remains so in the 0(3) non-linear a-model even on 
coarse lattices. We simulated the FP-action and the standard action at correlation 
length ~ 3 (/3pp=0.7 and /?5t=1-18, respectively). We have calculated the two- 
point function on a 24 x 24 periodic lattice. This lattice is large enough to avoid 
rotation symmetry breaking from the box itself (which is a physical infrared effect) 
up to distances ~ 9. In Figs. 8,9 the measured points are connected by piece-wise 
straight lines to lead the eyes. No symmetry breaking effects can be seen in the 
two-point function of the FP-action, while the lattice structure shows up clearly 
when the standard action is simulated. 

5 Surprisingly, even for an extremely coarse lattice, L/a = 3, (at (3fp = 0.98) one obtains 
g(2L) = 1.2626(11), still with no sign of lattice artefacts. 
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3 RG transformation at finite (3 



The FP-action is the perfect classical action, but we need the perfect quantum 
action. Since the FP-action seems to have small cut-off effects even at intermediate 
(3 values, it is a good starting point to find the RT. Assume that at some (3pp 
value the cut-off effects in the predictions of the FP-action (3fpA*{S) are negligible. 
Performing a RG step, the size of the cut-off effects for the blocked action will be the 
same (i.e. negligible), but the correlation length will be a factor of 2 smaller. This is 
the basic step we have to do in searching for the perfect action. In the following we 
shall illustrate that it is relatively easy to determine and parametrize the blocked 
action, and this blocked action remains short ranged if the RG transformation is 
chosen properly. 



3.1 The RG transformation 

We shall perform one RG step starting with a FP-action (in its parametrical form 
as given in Table 4) at (3fp = 1.0. At this point the correlation length is O(10). We 
use the block transformation defined in eq. @. There is no reason to assume that 
the optimal value of the parameter P is /3fp ■ 2 as it is in the case for large @fp- 
Actually, we expect that P op t goes to a constant rather than to zero for small (3fp- 
Preliminary runs indicated that k = 2.5 is close to the optimal choice, so we took 
P = f3pp ■ 2.5. In this case the range of interaction of the blocked action turned out 
to be essentially the same as that of the FP-action. For this reason it is sufficient 
to start with a 10 x 10 lattice and block it down to 5 x 5. At this small lattice sizes 



the Ferrenberg-Swendsen [20| technique becomes especially effective and we used it 
extensively. 



3.2 Determination of the properties of the blocked action 



In a MC calculation it is relatively easy to determine the change of the action under 
the change of the configuration. This information will be enough for us. Let us 
introduce the notation 



T(R,S) = £ 
Eq. (Rf) gives then 



PR nB £ S n -\nY N (p\ J2 S r 



e -f}'[A'(R)-A'{Ho)} _ / e T(R,S)-T(R ,S) 



Ro 



where 



«?>Ro 



J e -i3 FP A*(S)+T(R ,S)Q 
J e -/3 FP A*(S)+T(Ra,S) 



(34) 

(35) 
(36) 
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r 


Mr) 


r 


Mr) 


(1,0) 


0.549(1) 


(2,0) 


0.0037(14) 


(1,1) 


0.167(1) 


(2,1) 


0.0027(13) 






(2,2) 


0.0002(13) 



Table 5: After one RG step on the FP-action at (3 = 1.0 the new action f3'A! 
is parametrized according to eq. (12'). The leading couplings in f3'p'(r) are given 
above. The conventions in eq. ( pEUD give /3' = 0.92(2). 

and {R} and {Ro} are two different block-configurations. 

The difference T(R, S) — T(R , S) depends on the block-averages 

fn B = £ S n . (37) 

n&iB 

Simulating the effective action defined in eq. (|36|) with some fixed {Ro} and 
storing {f„ s } after every sweep, the expectation value in eq. (|35| ) can be calculated 
for any configuration {R}. In practice {R} should be close to {Ro}, otherwise the 
statistical errors can not be controlled. 

Specifically, one can chose {Ro} to be the trivial configuration (all the spins point 
in the 1 st direction). In this case ^4'(Ro) = (by definition; the constant generated 
by the RG transformation does not interest us), and we can calculate ^4'(R) for 
configurations which are close to the trivial one. By taking a parametrization for 
A! in the form of eq. (12') we get 

P'p'(r) = P 2 Q L (r), r^0, (38) 

where (N = 3) 
The conventions 

£//(r) = 0, Ep'(^ 2 = - 4 (40) 

r r 

determine then p'(0) and /3', respectively. We summarized the results in Table 5. 
It is interesting to observe that the structure of p'{r) is quite similar to p(r). Most 
importantly, it is similarly short ranged. 

We can now take configurations {R}, where a single spin is rotated by an angle 
in a trivial background. The corresponding value of (3' A' is given in the second 
column of Table 6. As in the case of the FP, A*, one observes a simple power-like 



20 





one spin 


rotated 


two spins rotated 


d 


P'A 


1.432 d 2 


P'A 


direct (2,0) interaction 


0.09967 


0.01427(1) 


0.01423 


0.02850(2) 


0.000037(28) 


tt/4 


0.8780(7) 


0.8833 


1.753(2) 


0.0029(23) 


3tt/8 


1.953(3) 


1.987 


3.91(1) 


0.000(14) 


tt/2 


3.42(1) 


3.53 


6.91(6) 


-0.08(6) 


5tt/8 


5.22(3) 


5.52 






3tt/4 


7.26(6) 


7.95 






2.8 


9.31(9) 


11.23 






3.0 


9.8(1) 


12.89 







Table 6: One and two spins (at a separation (2,0)) are rotated in a trivial back- 
ground. The two-spin result (4 th column) is almost exactly twice of the one spin 
result, indicating that the direct (2,0) interaction is small. The two leading couplings 
of p', giving 1.432 dominate the action even at large angles. 

behaviour in which justifies a parametrization in the form of eq. (12'). Actually, 
the two important couplings in p' which give a contribution 1.432 d 2 (third column 
in Table 6) dominate the results even at large angles When two spins, separated 
by the vector (2, 0), are rotated together, the measured action (column 4 in Table 6) 
is very closely twice of the result for a single rotated spin. The direct (2, 0) coupling 
is small. 

Since A' is short ranged, for the general parametrization we might consider 
the type of couplings used for the FP-action before. Actually, as we shall see, 
the FP-action (as parametrized in Table 4.) f3' FP A* describes P'A quite well with 
P'pp ~ 0.85. In applications the blocked action P'A will be used in a MC simulation, 
therefore we want a parametrization which gives the value of the action precisely 
for typical configurations. We generated 300 configurations by MC using the FP- 
action with P' FP = 0.85. We calculated the value of the blocked action P'A on 
these configurations using the technique discussed before. The type of the couplings 
kept for the fit where those used for the FP-action, plus the (2, 0), (2, 1) and (2, 2) 
two-spin couplings and their powers up to i? 6 . In order to see the quality of the fit 
and to separate in the fitting error the statistical (coming from MC determination 
of the action) and systematical (coming from parametrization) errors we kept those 
configurations only where the statistical error was below a certain cut value. Fig. 10 
shows the error of the fit against the value of the action. The average value of the fit 
error is 0.004. If the (2, 0), (2, 1) and (2, 2) couplings are left out, the average fit error 
is increased to 0.010. If needed, the quality of the fit can be increased by including 
further couplings (and increasing the statistical precision). In our feasibility study 
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we did not investigate this problem further. We want to mention two additional 
points only. The FP-action (3' FP A* with f3' FP = 0.848 gives a quite reasonable fit 
also with an average fit-error of 0.022. The second remark concerns the value k = 2.5 
which makes the blocked action somewhat more local and closer to the FP-action 
than K = 2we used at (3 — oo. At K — 2, for example, /3'p'(2, 0) = 0.011(1), a factor 
of ~ 3 larger than at k = 2.5. Similarly, the fit-errors (analogous to those discussed 
above) are increased by a factor of ~ 2 if k = 2 is taken. 

Acknowledgement. The authors are indebted for the useful discussions to M. Blatter, 
R. Burkhalter, G. Mack, K.-H. Mutter, A. Pordt, K. Schilling, P. Weisz, U. Wolff 
and U. Wiese. We thank M. Liischer for raising our attention to a minor flaw in the 
proof of the statement on classical solutions in the preprint version of this paper. 



A The Lambda parameter 

The A parameter for some specific lattice regularizations (in units of A$t corre- 
sponding to the standard action, for example) has been obtained earlier ET 



We give here the result for the generic case, for an action of the form given by 



eq. (P). 



The simplest way to calculate the ratio of the A parameters is to introduce a 
chemical potential h [^5] coupled to a Noether charge of the O(A) symmetry and 
use the fact that the free energy f{h) depends on h in a universal way ||25|| . 

Technically, the chemical potential is a constant, imaginary gauge potential: 
— > ih5^ q, where q is a generator of an O(A) rotation. We choose a rotation in 
the 1,2 plane: 



/ -i 
i 







o\ 




0/ 



(41) 



In the presence of the chemical potential the only change in the action is in the 
scalar product of two spins: 



Sr, S 



>->n" O r 



SriSr 



n J n+r ' J n.+r ~ J n J n+r (42) 

+ ihrQ ^S^ ^S^ r — ^S^ r ^ + — h r Q (s^Sn+r + S^S^+, r ^ + . . . 

(In fact, h is here ah and the omitted terms do not contribute in the continuum 
limit.) 



Using the representation in eq. (0) and expanding in n, after straightforward 
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calculations we obtain the following expression for the free energy: 

h 2 1 r ( 1 

f( h ) = ~7T + 7) / ln { P(q) + h2 + - Yl P( r ) r o cos qr-h 2 ^2 c(n 1} n 2 ; n 3 , n 4 ) 

Z 9 1 J <1 { A r m,...,n 4 



8r r' Q cos qd sin — sin — + Tq(1 — cosgr') + r' (1 — cosgr) 



(43) 



1 



+-(JV - 2) \n\ p(q) + h 2 -h 2 J2 c (- • •) r o(! ~ cosgr') + r' \l - cosgr) 

2 h { m,...,rH 



Expanding in h and using the fact that f(h) — /(0) is the same for any two 
regularizations, one obtains: 



where 



- + Q = -- + Q' 

9 9 



(44) 



Q= /^{iV-2-iEp(^ 2 (l-cosgr) 

a / \ / /s . . qr . q'r 

—4 c(ni, n.3, 7i4j(r • r J cosga sm — sin — 
ni,...,n 4 2 2 

— ~{N — 1) c(ni, n 2 ; n 3 , n 4 ) r 2 (l — cos gr') + r' 2 (l — cos gr) 

ni,...,n4 



(45) 



Here we introduced the notations for the relative coordinates in c(ni, n 2 ; n 3 , n 4 ): 



and used the 90° symmetry to replace e.g. 2ror[) by the scalar product (r • r'). 
Eq. (0) gives for the ratio of two A parameters: 



A' 



A =6XP l N-2 



2tt 



(Of - Q) 



(47) 



Note that only the 2- and 4-spin interactions contribute when the action is 
written in the form of eq. fll~2|) . 



In ref. |22| the A parameter has been calculated for a particular form of the 
4-spin interactions. Unfortunately we do not agree with the result quoted there. 

Since in numerical simulations one has to work with a restricted set of couplings, 
it is more informative to have the A parameter for that particular set. For the 
couplings used in our simulations the ratio of this A parameter to that of the standard 
lattice action is 8.17, as cited in eq. ([32]). 
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The A-parameter is rather sensitive to the coefficients in the action. We did not 
attempt to calculate it for the whole set of couplings in *4*(S), although this would 
be a feasible task. It is remarkable that while A$t is anomalously small compared 
to those in a continuum regularization (e.g. Aj^/Ast = 27.21), App is much closer 
to its continuum counterparts. 
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Figure captions 



Figure 1. The cut-off dependence of m(2L)2L for fixed value of m(L)L = 1.0595, 
taken from |2[. The values of L/a are indicated in the plot. The curve is the result of 
a fit with a second order polynomial in (a/L) 2 . The open box shows the extrapolated 
value. 

Figure 2. Flow of the couplings under RG transformation in the O(N) non-linear 
(j-model. 

Figure 3. The same as Fig. 1 with addition of two points generated by using the 
FP-action. These data show no cut-off dependence and agree with the value ex- 
trapolated to a/L = for the standard action. 

Figure 4. Dependence of the two-spin interaction coefficient p(r) from the distance 
r = (r , 0) for different choices of the block spin transformation. The circles and 
boxes correspond to the optimal choice with k = 2, and k = oo, respectively. 

Figure 5. The energy spectrum of the lattice regularized free field theory for different 
lattice actions: the FP (exact) (solid line), the approximation where the two leading 
couplings of the FP are kept only (dotted line), the standard action(dashed-dotted 
line) and the Symanzik tree level improved action (dashed line). Note that in the 
last case the eigenvalue becomes complex above \qi\ ~ 1.85. For these values the 
real part is plotted. 

Figure 6. The value of A*(S) when one of the spins is rotated by an angle $ with 
respect to a trivial background of parallel spins. The curve is the contribution from 
the quadratic and quartic terms in eq. (12'), Ai(&) = 3.2403($ 2 /2) - 0.0476(tf 2 /2) 2 . 

Figure 7. The quality of the parametrization of A* in Table 1 is shown in this scatter 
plot. The deviation from the true value of the action is plotted against the value of 
the action itself for « 500 configurations. 

Figure 8. The spin-spin correlation function as measured with the standard action 
at (3st = 1-18 corresponding to a correlation length ~ 3. The points are connected 
by piece-wise straight lines. 

Figure 9. The same as in Fig. 8 for the FP-action at (3fp = 0.7 corresponding to 
the same correlation length. 

Figure 10. The quality of the parametrization for the action (3' A' obtained after a 
block transformation on the FP-action at f3pp = 1.0 is shown here. The fit included 
the interactions displayed in Table 4 plus the (2,0), (2,1) and (2,2) couplings. 
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